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THE ENTROPY METHOD FOR REACTION-DIFFUSION SYSTEMS 
WITHOUT DETAILED BALANCE: 

FIRST ORDER CHEMICAL REACTION NETWORKS 

KLEMENS FELLNER, WOLFGANG PRAGER, BAO Q. TANG 


Abstract. In this paper, the applicability of the entropy method for the trend towards equilibrium for 
reaction-diffusion systems arising from first order chemical reaction networks is studied. In particular, 
we present a suitable entropy structure for weakly reversible reaction networks without detail balance 
condition. 

We show by deriving an entropy-entropy dissipation estimate that for any weakly reversible network 
each solution trajectory converges exponentially fast to the unique positive equilibrium with computable 
rates. This convergence is shown to be true even in cases when the diffusion coefficients of all but one 
species are zero. 

For non-weakly reversible networks consisting of source, transmission and target components, it is 
shown that species belonging to a source or transmission component decay to zero exponentially fast 
while species belonging to a target component converge to the corresponding positive equilibria, which 
are determined by the dynamics of the target component and the mass injected from other components. 
The results of this work, in some sense, complete the picture of trend to equilibrium for first order 
chemical reaction networks. 


1. Introduction and Main results 

This paper investigates the applicability of the entropy method and proves the convergence to equi¬ 
librium for reaction-diffusion systems, which do not satisfy a detailed balance condition. 

The mathematical theory of (spatially homogeneous) chemical reaction networks goes back to the 
pioneer works of e.g. Horn, Jackson, Feinberg and the Volperts, see [Fei79, Fei87, FH, Hor72, Hor74, 
HJ72, Vol, VVV] and the references therein. The aim is to study the dynamical system behaviour of 
reaction networks independently of the values of the reaction rates. It is conjectured since the early of 
1970s that in a complex balanced system, the trajectories of the corresponding dynamical system always 
converge to a positive equilibrium. This conjecture was given the name Global Attractor Conjecture by 
Craciun et al. [CDSS]. The conjecture in its full generality is - up to our knowledge - still unsolved so 
far, despite many attempts have been made by mathematicians to attack this problem. 

From the many previous works concerning the large time behaviour of chemical reaction networks, the 
majority of the existing results considers the spatially homogeneous ODE setting. The PDE setting in 
terms of reaction-diffusion systems is less studied. Also detailed quantitative statements like, e.g. rates 
of convergence to equilibrium, constitute frequently open questions even in the ODE setting. 

Our general aim is to prove quantitative results on the large-time behaviour of chemical reaction 
networks modelled by reaction-diffusion systems. In the present work, we study reaction-diffusion systems 
arising from first order chemical reaction networks and show that all solution trajectories converge 
exponentially to corresponding equilibria with explicitly computable rates. 

Our approach applies the so called entropy method. Going back to ideas of Boltzmann and Grad, 
the fundamental idea of the entropy method is to quantify the monotone decay of a suitable entropy 
(e.g. a convex Lyapunov) functional in terms of a functional inequality connecting the time-derivative 
of the entropy, the so called entropy dissipation functional, back to the entropy functional itself, i.e. to 
derive a so called entropy entropy-dissipation (EED) inequality. Such an EED inequality can only hold 
provided that all conservation laws are taken into account. After having established an EED inequality 
and applying it to global solutions of a dissipative evolutionary problem, a direct Gronwall argument 
implies convergence to equilibrium in relative entropy with rates and constants, which can be made 
explicit. 
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By being based on functional inequalities (rather than on direct estimates on the solutions), a major 
advantage of the entropy method is its robustness with respect to model variations and generalisations. 
Moreover, the entropy method is per se a fully nonlinear approach. 

The fundamental idea of the entropy method originates from the pioneer works of kinetic theory and 
from names like Boltzmann and Grad in order to investigate the trend to equilibrium of e.g. models of 
gas kinetics. 

A systematic effort in developing the entropy method for dissipative evolution equations started not 
until much later, see e.g. the seminal works [Tos, TV, CJMTU, AMTU, DVOl] and the references therein 
for scalar (nonlinear) diffusion or Fokker-Planck equations, and in particular the paper of Desvillettes and 
Villani concerning the trend to equilibrium for the spatial inhomogeneous Boltzmann equation [DV05]. 
The derivation of BED inequalities for scalar evolution equations is typically based on the Barky-Emery 
strategy (see e.g. [CJMTU, AMTU]), which seems to fail (or be too involved) to apply to systems. 

The great challenge of the entropy method for systems is, therefore, to be able to derive an entropy 
entropy-dissipation inequality, which summarises (in the sense of measuring with a convex entropy func¬ 
tional) the entire dissipative behaviour of solutions to a (possibly nonlinear) dynamical system to which 
the EED inequality shall be applied to. Preliminary results based on a (non-explicit) compactness- 
contradiction argument in 2D were obtained e.g. in [Gro, GGH, GH] in the context semiconductor 
drift-diffusion models. 

The first proof of an EED inequality with explicitly computable constants and rates for specihc nonlin¬ 
ear reaction-diffusion systems was shown in [DE06] and followed by e.g. [DE07, DF08, DF15, FLT, MHMj. 
The application of these EED inequalities to global solutions of the corresponding reaction-diffusion 
systems proves (together with Csiszar-Kullback-Pinsker type inequalities) the explicit convergence to 
equilibrium for these reaction-diffusion systems. 

We emphasise that all these previous results on entropy methods for systems assumed a detailed 
balance condition and, thus, features the free energy functional as a natural convex Lyapunov functional. 

A main novelty of the paper lies in demonstrating how the entropy method can be generalised to first 
order reaction networks without detailed balance equilibria. In particular we shall consider firstly weakly 
reversible networks and secondly even more general composite systems consisting of source, transmission 
and target components (see below for the precise definitions). 

We feel that it is important to point out that while there are certainly many classical approaches by 
which linear reaction-diffusion systems can be successfully dealt with, our task at hand is to clarify the 
entropic structure and the applicability of the entropy method for linear reaction networks as a first step 
before being able to turn to nonlinear problems in the future. See [DFT] for such a generalisation of the 
method to nonlinear reaction-diffusion systems satisfying the so-called complex balance condition (see 
Definition 1.3 below). 

The goal of this present work is to prove the explicit convergence to equilibrium for the complex 
balanced and more general reaction-diffusion systems corresponding to first order reaction networks. To 
be more precise, we study first order reaction networks of the form 

i^j = l,2,...,N, (AT) 

aji 


Eigure 1 . A first-order chemical reaction network 

where Si,i = 1, 2,..., V, are different chemical substances (or species) and aij, Oji > 0 are reaction rate 
constants. In particular, Oij denotes the reaction rates from the species Sj to Si. 

Eirst order reaction networks appear in many classical models, see e.g. [Smo, Rot]. More recently, 
first order catalytic reactions are used to model transcription and translation of genes in [TVO]. The 
evolution of the surface morphology during epitaxial growth involves the nucleation and growth of atomic 
islands, and these processes may be described by first order adsorption and desorption reactions coupled 
with diffusion along the surface. A first order reaction network can also be used to describe the reversible 
transitions between various conformational states of proteins (see e.g. [MGetal]). RNA also exists in 
several conformations, and the transitions between various folding states follow hrst order kinetics (see 
[BRetal]). 

In the present paper, we investigate the entropy method and the trend to equilibrium of reaction- 
diffusion systems modelling first order reaction networks with mass action kinetics. More precisely, we 
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shall consider the reaction network M in the context of reaction-diffusion equations and assume that 
for all i the substances Si are described by spatial-temporal concentrations Ui{x,t) at 

position X G fl and time t >0. Here, H shall denote a bounded domain H C K" with sufficiently smooth 
boundary dfl (that is d^l G to avoid all difficulties with boundary regularity, although the below 

methods should equally work under weaker assumptions) and the outer unit normal v{x) for all x G dVl. 
Due to the rescaling x |r2|^/"a:, we can moreover consider (without loss of generality) domains with 
normalised volume, i.e. 

|D| = 1. 

In addition, we assume that each substance Si diffuses with a diffusion rate > 0 for all i = 
1,2,..., N. Finally, we shall assume mass action law kinetics as model for the reaction rates, which leads 
to the following linear reaction-diffusion system: 

{ Xt = DAX + AX, xGn, t>0, 

dyX = 0, X G dfl, t > 0, (1.1) 

X(a:,0) = Xo(x) > 0, x G n. 


where X{x,t) = \ui{x,t),U2{x,t),... ,UN{x,ty\^ denotes the vector of concentrations subject to non¬ 
negative initial conditions Xq{x) = [ui_o(a:) > 0, U2,q{x) > 0,..., un,o{x) > 0]^, D — diag(di, d2, ■ ■ ■, (In) 
denotes the diagonal diffusion matrix and the reaction matrix A = (aij) G satisfies the following 

conditions: 


aij>0, foralH^j, ■i,j = l,2,...,X, 

for all j = 1, 2,..., iV. 


( 1 . 2 ) 


The conditions (1.2) on the reaction matrix A imply in particular that the vector (1,1,...,!)^ con¬ 
stitutes a left-eigenvector corresponding to the eigenvalue zero. Together with homogeneous Neumann 
boundary conditions this implies that solutions to (1.1) admit the following conservation of total mass: 

N 

Ui{x,t)dx = 'S^ / Uifi{x)dx =: M > 0, for alH > 0, (1.3) 

i=i 



where M > 0 is the initial total mass, which we shall assume positive. 

If X{x,t) = X{t), then system (1.1) reduces to the corresponding space-homogeneous ODE model. 
Independently of PDE- or ODE-setting, we recall the following definitions of equilibria from e.g. [HJ72, 
Fei79, VVV]. 


Definition 1.1 (Homogeneous Equilibrium). 

A state Xoo = (ui^oo, M 2 ,oo, ... ,un,oo) is called a homogeneous equilibrium or shortly equilibrium of the 
first order reaction network Af if AXoo = 0. 


Definition 1.2 (Detailed Balance Equilibrium). 

A positive equilibrium state Xoo = (mi,oo, R 2 ,oo,..., my.oo) > 0 is called a detailed balance equilibrium 
for the reaction network Af if a positive reaction rate constant atj > 0 for i ^ j implies also a positive 
reversed reaction rate constant Oji > 0 and that the forward and backward reaction rates balance at 
equilibrium, i.e. 

The reaction network Af is called to satisfy the detailed balance condition if it admits a detailed balance 
equilibrium. 


Definition 1.3 (Complex Balance Equilibrium). 

A positive equilibrium state Xoo = (ui.oo'«2,oo5 ■. ■, waTjOo) > 0 zs called a complex balance equilibrium 
for the reaction network Af if for all k = 1,2,..., N, the total in-flow into the substance Sk balances in 
equilibrium the total out-flow from Sk to all other substances Si, i.e. 

^ ^ e^ki^i.oo — ( ^ ^ e^jk j'^kjOO- 

{l<i<N-. aM>0} ajk>0} ' 

The reaction network Af is called complex balanced if it admits a complex balance equilibrium. Moreover 
for complex balanced chemical reaction networks, all equilibria are complex balanced, see e.g. [Hor72]. 

Example 1.1 (Detailed balance equilibria are complex balance equilibria). 

It is easy to see that detailed balance equilibria are also complex balance equilibria while the reverse does 
not hold in general, even for reversible networks. For example, consider the reaction network in Figure 
2, where all reaction rates constants Oij > 0 are assumed positive and the network is thus fully reversible. 
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Figure 2. A reversible network 


The corresponding reaction-diffusion system with homogeneous Neumann boundary conditions 

dtui - diAui = -(021 + 031)01 + 012M2 + 01303, 
dtU2 - d2Au2 = O21O1 - (012 + 032)02 + O23O3, 
dtU3 - dsAua = 03101 + 032O2 - (013 + 023)03, 
d„ui = di,U2 = dvU3 = 0. 

exhibits the constant equilibrium Xoo = (oi_oo, 02 , 00 , 03 , 00 ) satisfying AXoo = 0 , i.e. 


{ ^12^2,00 “1“ ^ 13 ^ 3,00 — (021 + 031)01 ,00 j 

021^1,00 H” ^23^3,00 = (012 + 032)02 ,00 j (1.5) 

^31^1,00 H” ^32^2,00 = (oi3 + 023)03 

,CXD ; 


which has a unique nontrivial solution once the mass conservation (1.3) is taken into account. 

According to Definition 1.3, it is clear that system (1.5) constitutes a complex balance equilibrium 
for all reaction rate constants Oij > 0. For Xca to be a detailed balance equilibrium, however, it is 
additionally necessary that 

{ 01202,00 = 02101^00, 

02303,00 = 03202 , 00 , ( 1 - 6 ) 

^31^1,00 — ^ 13 ^ 3,005 

which obviously implies (1.5). Yet the equations (1.6) can only have a solution if 


012 • 023 • 031 
O21 • O32 • Oi3 


(1.7) 


holds; in other words if the product of the reaction rate constants multiplied in the clockwise sense 
of the above reaction network graph equals the product of the reaction rate constants multiplied in the 
counterclockwise sense. The condition (1.7) is thus necessary and sufficient for system (1.4) to admit a 
detailed balance equilibrium. 


Remark 1.1 (General definition of detailed and complex balance). 

The concepts of detailed balance and complex balance are also defined for general higher order chemical 
reaction networks, see e.g. [Hor72]. For simplicity, we stated here the definition corresponding to the 
first order network A/". In general, one can roughly say that a state X^o is called a complex balanced 
equilibrium if at equilibrium the total in-flow to each specie Si is equal to the total out-flow from Si. 

Remark 1.2 (Detailed balance and reversibility). 

It follows from Definition (1.2) that if M satisfies the detailed balance condition, then it is also reversible 
in the sense that for any reaction Si —>■ Sj also the reverse reaction Sj —>■ Si takes place. 


The set of complex balanced systems is much larger than the one of detailed balance systems. Horn 
already gave necessary and sufficient conditions for a network to satisfy the complex balance condition 
in [Hor72]. For convenience of the reader, we present in the following the associated definitions of 
directed graphs as representations of reaction networks. The image of the associated graphs will also 
help following some of our main estimates. 

A directed graph G corresponding to a given reaction network M is defined by considering the sub¬ 
stances Si,i = 1,2,... ,N, as the N nodes of G, which are connected for alH j = 1, 2,..., by an 

dji 

edge with starting node Si and finishing node Sj if and only if the reaction Si —> Sj occurs with a 
positive reaction rate constant aji > 0. 
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Definition 1.4 (Linkage classes partition of a first order reaction network, Connected networks). 

A linkage class C of a first order network M is a maximal set of connected substances, i.e. Si, Sj G C 
implies that Si and Sj are connected (in the sense that there exist Si = S ^, Sr^ ■ ■ ■, Sr^-i, Sr^ = Sj such 
that for each 1 < £ < k — 1, either the reaction S^ —>■ or Sr/.^^ —>■ happens) but Si G C and 

Sj ^ C implies that Si and Sj are not connected. 

If a reaction network consists only of one linkage class, we shall call such a network connected. 

Definition 1.5 (Weak reversibility of a first order reaction network). 

A first order reaction network Af is called weakly reversible if for any reaction Si —>■ Sj with i j, there 
exists a chain of reactions Sj = Sj., —>■ Sj.^ Sj,. = Si where Sj,, Sj .^, ■ • ■, Sj, are other chemical 

substances of Af . 

If a reaction network Af is weakly reversible, then we also call the corresponding directed graph G 
weakly reversible. 

Definition 1.6 (Strongly connected components of a directed graph). 

A subgraph H C G of a directed graph G is called a strongly connected component if for any two nodes 
Si,Sj in H, we can find a path from Si to Sj of the form Si ^ Si, ^ ^ Si, —>■ Sj with all 

Si,, Si 2 ,..., Si, belonging to H. 

We call a first order reaction network Af strongly connected when its corresponding graph G is strongly 
connected. 

Remark 1.3 (Partition of weakly reversible first order reaction networks Af into disjoint strongly con¬ 
nected components/subnetworks). 

Firstly, it follows directly from Definition 1.4 that any first order reaction network Af can be uniquely 
partitioned into a pairwise disjoint union of linkage classes and each linkage class C constitutes a con¬ 
nected subnetwork Afc ■ In particular, for a weakly reversible first order reaction network Af, each linkage 
class C forms a connected weakly reversible subnetwork Afc and it is straightforward to show that the 
directed graph corresponding to Afc is strongly connected according to Definition 1.6. (Consider that 
for all reactions being part of the connection between Si,Sj G Afc, the weak reversibility implies the 
existence of a returning chain of reactions. Thus, there exist chains of reactions connecting Si to Sj 
and vice versa.) black Secondly, any directed graph G can be partitioned into a pairwise disjoint union 
of strongly connected components, all of which are weakly reversible according to Definition 1.5. Note 
that these strongly connected components can still be connected via “non-weakly-blackreversible” reactions 
(see e.g. Figure 3). Therefore, for general directed graphs, multiple strongly connected components may 
constitute one linkage class. However, if the directed graph G is additionally weakly reversible, then each 
strongly connected component has to constitute exactly one linkage class since otherwise we have already 
seen that weakly reversible subnetworks Afc corresponding to one linkage class C are strongly connected. 

Thus, for weakly reversible first order reaction networks Af, the partition of linkage classes is identical 
to the partition of strongly connected components of the corresponding directed graphs. 

Therefore, with a marginal abuse of notation, we will use the terminology “strongly connected compo¬ 
nent” or “strongly connected subnetwork” both for such a connected weakly reversible first order reaction 
subnetwork Afc o-nd its corresponding strongly connected subgraph/component. 

Remark 1.4 (Linkage classes of first order reaction networks can be treated independently). 

For first order reaction networks, each node represents exactly one substance. Thus, any linkage class of 
a first order reaction network can be treated independently from the others. In particular, all the strongly 
connected components of a weakly reversible first order reaction network can he treated independently 
since these subnetworks form different linkage classes. 

For higher order reaction networks, where the nodes of the corresponding graphs are so-called com¬ 
plexes consisting of multiple substances, this is not necessarily true since one substance might need to be 
represented by different nodes. 

Because of Remarks 1.3 and 1.4, we will consider in Section 2 weakly reversible first order networks 
partitioned into strongly connected first order reaction subnetworks Afc, and each strongly connected 
component Afc can (w.l.o.g) be treated independently. In Section 3, we will consider (w.l.o.g) connected 
reaction networks Af consisting of one linkage class, yet we shall not assume weak reversibility. Hence 
the corresponding directed graphs are not strongly connected and may consists of multiple strongly 
connected components, but the underlying undirected graphs are connected (see e.g. Figure 3). 

Lemma 1.1 (Strongly connected networks, irreducible reaction matrices and complex balance equilib¬ 
ria). 

For any first order reaction network Af the following statements are equivalent: 
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• The first order reaction network Af is strongly connected. 

• The corresponding reaction matrix A of Af is irreducible. 

• The first order reaction network Af is complex balanced and for any positive mass M > 0 (as set 
by the conservation law (1-3) and there exists of a unique, positive complex balance equilibrium 
Xoo = (ui.oo, U2,oo, • ■ •, un,oo) > 0 o/ System (1.1), which satisfies 


AXoo = 0, 

Efci = M > 0. 


( 1 . 8 ) 


Proof. The equivalence of strong connectivity for first order networks and irreducibility of the reac¬ 
tion matrix A follows e.g. from [SenSl, Definition 2.1, page 46] and [Min88, Theorem 3.2, page 78]. 
Next, the Perron-Frobenius theorem implies for any irreducible reaction matrix A and any positive mass 
'u-i.co = M > 0 the existence of a unique positive equilibrium, see e.g. [Sen81, Per07] and Lemma 
2.2 below. This equilibrium satisfies AX^c — 0 and is thus a complex balance equilibrium according 
to Definition 1.3. Hence, the strongly connected first order reaction network Af is complex balanced 
(independently of the value of M). Finally, Lemma 2.2 below implies that strongly connected first order 
reaction networks possessing unique positive equilibrium (for fixed M > 0) have irreducible reaction 
matrices A. □ 


Remark 1.5 (Complex balanced higher order systems are necessarily weakly reversible). 

For higher order reaction network, it holds only true that systems with complex balance equilibrium are 
necessarily weakly reversible. Thus, weakly reversible systems constitute the more general class of reaction 
networks. 


Remark 1.6. The equilibrium Xoo in (1.8) is spatially homogeneous. Thus, it coincides with the equilib¬ 
rium for the corresponding spatially homogeneous ODE system Xt = AX of the reaction network given 
in Figure 1. In [And] or [SiMa], the authors proved that X{t) —> X^o as t —> -l-oo. However, the 
method used in this paper cannot be directly applied to prove the convergence to equilibrium for PDE 
system (1.1). 


The first main result of this paper concerns the convergence to equilibrium for weakly reversible 
reaction networks of the form displayed in Figure 1. Our method of proof applies the entropy method 
to prove explicit exponential convergence of solutions of system (1.1) to the unique equilibrium. 

As mentioned above, all previous results of explicit FED inequalities (see e.g. [DF06, DF07, DF08, 
DF15, FLT, MHM]) considered reaction-diffusion systems satisfying a detailed balance condition. 


In the current paper, we shall show that the following quadratic relative entropy between any two 
solutions X = {ui,... ,un) and Y = (vi,..., vn) 

£(X\Y){t)=Y, [ ^-^dx (1.9) 

is an entropy functional, see Lemma 2.3 below, which is the first key result of this paper. 

In particular, we can consider the special case Y = X^o for such an entropy functional. By using the 
linearity of first order systems, it is then straightforward to check (by using (1.2) and AX^o = 0) that 
the quadratic relative entropy towards an equilibrium state Xao, i.e. 

£{X - XooIXoo) = V / (110) 

i=l 

is equally an entropy functional, which decays monotone in time according to the following explicit form 
of the entropy dissipation functional — Xod^oo) = —'D{X — Xoc\Xoc)‘- 


N 


V{X-Xoo\Xoo) = 2 Y,d^ 


\X{lli 'Ui^oo') I 


dx 


N 


N ^ 1^ 19 N 


f ( 

( Ui Ui^QQ 

^j ,cso N 

/o' 

^ ^2,00 

J 


dx > 0 


(1.11) 


, f . f f u 

— 2 y ^ di I dx ^ H” ^ij^j,oo / / ( 

i=i dn dn \Ui, 


Ui Ui 


oo '^J,oo 


dx 


= V{X\Xoo) = --£{X\X^) > 0. 
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The dissipative structure of the quadratic relative entropy towards equilibrium (1.10) is a special 
cases of generalised relative entropies discussed e.g. in [Per07, Chapter 6]. The entropy functional (1.9), 
i.e. the observation of the dissipativeness of the relative entropy between any two solutions, is however 
related to a general property of linear Markow processes, which was recently shown in [FJ16]. 

With the help of the explicit form of entropy dissipation (1.11), we are able to show (in Lemma 2.4 
below) an entropy-entropy dissipation inequality of the form 

V{X-X^\X^)>\S{X-X^\X^), (1.12) 

where A > 0 is an explicitly computable constant. Once the EED inequality (1.12) is proven, the 
statement of the first main theorem follows from a standard Gronwall argument, see Section 2 below: 

Theorem 1.2 (Exponential equilibration of weakly reversible first order reaction networks). 

Given a weakly reversible first order reaction network partitioned into linkage classes. Consider (w.l.o.g.) 
any corresponding strongly connected subnetwork Mc ■ Assume for Me that the diffusion coefficients di 
are positive for all i = 1,2,N, and the initial mass M is positive. 

Then, the unigue global solution to initial-boundary problem (1.1) converges exponentially to the unique 
positive equilibrium X^o = (ui^oo, ^ 2 , 00 , •. ■ ,un,oo), be. 

y [ ^ ^-xt y f 

^ Jn Wi.oo y Jn 

where the constant A > 0 depends explicitly on the reaction matrix A, the domain Q, the diffusion matrix 
D and the initial mass M. 

Remark 1.7 (Lyapunov functionals for ODE systems). 

For ODE systems, Lyapunov functionals have been mainly considered in the analysis of nonlinear ODE 
systems. Moreover, for nonlinear ODE systems, L^-type Lyapunov functionals are most commonly used 
in the study of the large-time-behaviour. Eor reaction-diffusion systems, however, L^-functionals are not 
useful for the entropy method and proving explicit convergence to equilibrium, since they do not measure 
the spatial diffusion in an exploitable way. 

We also remark, that while logarithmic relative entropy functionals of the form 

N 

Vx^{X){t) = ^ {ui{lnu^ - - 1) + Ui,oo) (1-13) 

were known to constitute a monotone decaying Lyapunov functional for complex balanced ODE reaction 
networks (see e.g. [HJ72, Gop, SiMa]^, up to our knowledge and somewhat surprisingly, no explicit 
expression of the entropy dissipation —dV/dt in complex balanced systems has been derived so far. 

We also refer the reader to e.g. [MiSi] for the stability of some mass action law reaction-diffusion 
systems, where the author used techniques ofui-limit sets along with the monotonicity of L^-type Lyapunov 
functional. 

Our results in this paper are significantly stronger in the sense that we show, by using the entropy 
method, the exponential convergence to equilibrium with computable rates. 

In addition and in comparison to co-limit techniques, the entropy method has also the major advantage 
of relying on functional inequalities rather than on specific estimates of solutions to a given system. 
Having such functional entropy entropy-dissipation inequalities once and for all established makes the 
entropy method robust with respect to model variations and generalisations. 

As example, it is the intrinsic robustness of the entropy method, which makes it possibly to also apply 
to non weakly reversible reaction networks, see Theorems l.f and 1.5 below. 

The assumption on the positivity of all diffusion coefficients in Theorem 1.2 is not necessary as such. 
As already shown in e.g. [DF07, FLT], the combined effect of diffusion of a specie and its weakly reversible 
reaction with other (possibly non-diffusive) species will lead to a indirect “diffusion-effect” on the latter 
specie. This indirect diffusion-effect can also be measured in terms of functional inequalities. Hence the 
exponential convergence to equilibrium still holds for systems with partial degenerate diffusion. 

Note that the indirect “diffusion transfer” and the convergence results of this paper resembles to 
some degree the framework of hypocoercivity for evolution equations like linear kinetic Fokker-Planck 
equations, see e.g. [Vil09, DMS, AAS]. However, while hypocoercivity typically requires the use of 
suitably constructed Lyapunov functionals, the indirect “diffusion-effect” can be entirely express in func¬ 
tional inequalities linking the relative entropy and the associated entropy dissipation functional. The 
entropy method present in this paper proves convergence to equilibrium essentially regardless of full- or 
degenerate diffusion matrices. 
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The exponential convergence for weakly reversible systems (1.1) with degenerate diffusion is stated in 
the following Theorem 1.3 to be proved in Section 2 below: 

Theorem 1.3 (Equilibration of linear networks with degenerate diffusion). 

Given a weakly reversible first order reaction network partitioned into linkage classes. Consider (w.l.o.g.) 
any corresponding strongly connected subnetwork Afc- Assume that the initial mass M is positive for 
Afc- Moreover, assume that at least one diffusion coefficient di is positive for some t = 1 , 2 ,..., N. 
Then, the solution to (1.1) converges exponentially fast to the unique positive equilibrium X^o = 

{ui 

,005*^2,00; • ■ • ; '^N ,Oo) • 

Jq '^1,00 Jq ^2,00 

with a computable rate X' > 0, which depends explicitly on A, 12, D and M. 

Remark 1.8 (Same results of linear ODE reaction networks). 

We remark that our approach can of course be adapted to equally apply to linear ODE reaction networks 
by eliminating the terms and calculations concerning .spatial diffusion. Thus, all the results of this paper 
hold equally for such linear ODE systems. 

As the second main result of this manuscript, we shall derive an entropy approach and prove conver¬ 
gence to equilibrium for reaction networks as in Eigure 1, for which the weak reversibility assumption 
does not hold. Eor first order reaction networks, this implies that the system is not complex balanced, 
or in other words, that equilibria are not necessarily positive. 

Due to the lack of positivity of equilibria, it follows immediately that the relative entropy used for 
weakly reversible systems is not directly applicable. In the following we proposed a modified entropy 
approach. At first, it is necessary to understand the structure of non weakly reversible reaction networks. 

We state here the necessary terminology and the main ideas. Since for any non weakly reversible 
linkage class, the associated directed graph G is connected (which means that the underlying undirected 
version of G is a connected graph) but not strongly connected, G consists of r > 2 strongly connected 
components, which we denote by Gi, G 2 ,..., G^. Then, we can construct a directed acyclic graph G^, 
i.e. G^ is a directed graph with no directed cycles as follows: 

a) G'^ has as nodes the r strongly connected components Gi, G 2 ,..., Gr, 

b) for two nodes Ci and Cj of G^, if there exists a reaction Ci 9 Sk Sf, S Cj with agk > 0, 
then there exists also the edge Ci —>■ Cj on G^. 

Due to the structure of G*^, its nodes, or equivalently the strongly connected components of G, can be 
labeled as one of the following three types: 

• A strongly connected component Ci is called a source component if there is no in-flow to Ci , i.e. 
there does not exist an edge Sk —>■ Sj where Sk ^ Ci and Sj € Ci. 

• A strongly connected component Ci is called a target component if there is no out-flow from Ci, 
i.e. there does not exist and edge Sk —>■ Sj where Sk G Ci and Sj ^ Ci. 

• If Ci is neither a source component nor a target component, then we call Ci a transmission 
component. 

Example 1.2. Consider the reaction network in Eigure 3. The depicted network has 4 strongly connected 

032 063 

02 -> >33 -*■ OQ 



Eigure 3. A non-weakly reversible reaction network consisting of four strongly con¬ 
nected components 

components Gi = { 5 'i,S' 2 }, G 2 = {•S'a}, G 3 = {S' 4 , 5 ' 5 }, G 4 = {.S'e}, where Gi is a source component, C 2 
is a transmission component and G 3 ,G 4 are target components. 
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By definition, each of the three types of strongly connected components is subject to a different 
dynamic, which can be written as follows: Let Ci be a strongly connected component and denote by Xi 
the concentrations within Ci. Moreover, denote by Ai the reaction matrix formed by all reactions within 
the component Ci. Then, we have 

• for a source component Ci\ 

dtX, - D,AX, = A,X, + Xr*, 

where summarises the out-flow from the source component Ci. 

• for a target component Ci'. 

dtX, - D,AX, = xr + AX^, 

where summarises the in-flow into the target component Ci. 

• for a transmission component Ci'. 

dtX, - D,AX, = Xr + A,X, + Xr\ 

where Xf^ , X°'^* are the in/out-flow of the transmission component Ci. 

In the dynamics of transmission and target components, the in-flow depends only on species which 
do not belong to Ci, so that X™ can be treated as an external source for the system for Ci. However, it 
may happen that Jy" contains inflow from species whose behaviour is not a-priori known. 

For acyclic graphs C^ , however, it is possible to avoid these difficulties, since the topological order of 
acyclic graphs allows to re-order the r strongly connected components Ci,C 2 , ... ,Cr in such a way that 
for every edge Ci —>■ Cj of it holds that i < j. This permits to study the dynamics of all components Ci 
sequentially according to the topological order and, when at times considering a transmission component 
(or later a target component) Ci, the required in-flow X™ contains only species whose behaviour is 
already known. 

Due to the structure of the network, it is expected that species belonging to source or transmission 
components are subsequently losing mass such that the concentrations decay to zero in the large-time 
behaviour as time goes to infinity. In contrast, the species belonging to a target component converge 
to an equilibrium state, which is determined by the reactions within this component and by the mass 
’’injected” from other components. 

Since the source and transmission components do not converge to positive equilibria, the relative en¬ 
tropy method used for weakly reversible networks is directly not applicable. Instead, for each component 
Ci we will modify the entropy method by introducing an artificial equilibrium state with normalised mass, 
which balances the reaction within Ci. The artificial equilibrium will allow us to consider a quadratic 
functional, which is similar to the relative entropy in weakly reversible networks and which can be proved 
to decay exponentially to zero. This result is stated in the following Theorem: 

Theorem 1.4 (Exponential decay to zero of source and transmission components). 

Given an arbitrary first order reaction network partitioned into linkage classes and consider (w.l.o.g.) 
any corresponding connected subnetwork Mc ■ Assume for Me that all diffusion coefficients di are positive. 

Then, for each Ci being a source or a transmission component of Me, there exist constants Ki > 0 
and Ai > 0 depending explicitly on Ai and fl such that, for any specie Si G Ci, the concentration ui of 
Si decays exponentially to zero, i.e. 

l|wf(*>-)lli 2 (Q) < Kie~^'\ for all t > 0. 

For a target component Ci, due to the in-flow Jq*”, the total mass of Ci is not conserved but increasing. 
Hence, Ci does not possesses an equilibrium as weakly reversible networks, which is explicitly given in 
terms of the reaction rates and the conserved initial total mass. 

However, since each target component is strongly connected and thus a weakly reversible reaction 
network with mass influx, there still exists a unique, positive equilibrium of Ci denoted by Xi^oo, which 
balances the reactions within Ci and has a total mass, which is the sum of the total initial mass of Ci and 
the total ’’injected mass” from the other components via the in-flow X™ (see Lemma 3.3). We emphasis 
that in general the injected mass is not given explicitly but depends on the time evolution of all the 
influencing species higher up with respect to the topological order of the graph G^. 

Since the equilibrium Xi^oo is positive, we can use again a relative entropy functional to prove the 
convergence of the species belonging to a target component to their corresponding equilibrium states. 

Theorem 1.5 (Exponential convergence for target components). 

Given an arbitrary first order reaction network partitioned into linkage classes and consider (w.l.o.g.) 
any corresponding conneeted subnetwork Me ■ Assume for Me that all diffusion coefficients di are positive. 
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Then, for all target components Ci = {Si^, Si^,..., Si,^ } of Me, where M is the number of species 
belonging to Ci , there exists a unique positive equilibrium state Xi^^o = (mii.oo, • ■ ■, miw., oo) (md the 
concentrations Ui^ of Si^ converges exponentially to the corresponding equilibrium value 

\\uii{t) - Ui,,oo||i 2 (o) < for all t > 0, 

with the constants Ki > 0 and Ai > 0 depending explicitly on Ai, and Di and on the equilibrium state 
Xi^oo • 


Remark 1.9. Note that by Lemma 3.3, the equilibrium state Xi^oo depends explicitly on the mass injected 
into the target component Ci, but that injected mass itself depends non-explicitly on the initial data and 
on the history of the reaction-diffusion network. 

Remark 1.10. We remark that in the same way as Theorem 1.3 generalises Theorem 1.2 to allow for 
degenerate diffusion matrices, it is equally possible to generalise Theorems l.f and 1.5 in the sense that it 
is sufficient to assume that for each target component there is at least one diffusion coefficient is positive. 
In particular, the proof of Theorem l.f holds independently from the entries of a non-negative diffusion 
matrices Di. 

Outline: The rest of the paper is organised as follows: In Section 2, we present the entropy method for 
weakly reversible networks and prove exponential convergence to the positive equilibrium. Non weakly 
reversible networks will be investigated in the Section 3. By using the structure of the underlying 
graphs, we are able completely resolve the large-time behaviour of all species belonging to such first 
order networks. 

We also remark that all constants in this manuscript are explicit in the sense that they are derived in 
constructive ways. However, since these constants are not optimal, we will denote them by using generic 
letters like Ki or Ai, etc. The issue of optimal rates and constants for the convergence is subtle, and can 
be investigated in future works. 

Notation: We shall use the shortcut / = /q f{x) dx, whenever |H| = 1, and || • || for the usual norm in 

i.e. 

Il/f = [ \fix)\^dx. 

Jn 

2. Strongly connected first order networks 

In this section, we consider strongly connected first order reaction networks M, for which the associated 
directed graph is strongly connected. This is w.l.o.g. by Remarks 1.3 and 1.4, since any weakly reversible 
first order reaction network can be partitioned into disjoint strongly connected components/subnetworks, 
which can be treated independently. 

Moreover, we recall that H C R" is a bounded domain with smooth boundary dll (say dU G C^+“) 
and normalised volume |r2| = 1 (w.l.o.g. by rescaling). Finally, we recall the system (1.1) 

{ dtX - DAX = AX, xGn, t>0, 
d,,X = 0, xGdfl, t>0, (2.1) 

A:(a:,0) = A:o(a:), x € Q, 

where X = \ui, U 2 ,. ■. ,un]’^ denotes the vector of concentrations, the vector Xq = [ui,o, U 2 ,o ,..., 
denotes the initial data, the diffusion matrix D = diag(di, (i 2 , • ■ •, ^at) and the reaction matrix A = 
(aij) G satisfies 

ioiaWiffj, i,j = l,2,...,N, 

[ajj = for all j = l,2,...,iV. 

Moreover, since Me is strongly connected, we know that the reaction matrix A is irreducible, see Lemma 
1.1. For the linear system (2.1), the existence of a global unique solution follows by standard arguments, 
see e.g. [Smo, Rot]: 

Theorem 2.1 (Global well-posedness of linear reaction-diffusion networks). 

For all given initial data Xq G , there exists a unique solution X G C([0, Tj; (L^(H))'^) n 




ENTROPY METHOD FOR REACTION-DIFFUSION SYSTEMS 


11 


L^{0,T; for all T > 0. Moreover, if Xq > 0 then X{t) > 0 for all t > 0. Finally, the 

solutions to (2.1) conserve the total mass (1.3) for all t > 0; 



N 

Ui{x, t)dx = / 


Uifi{x)dx =: M > 0, 


(2.3) 


where the initial mass M is assumed positive. 


Lemma 1.1 stated the equivalence to weak reversibility first order reaction networks and irreducibility 
of the reaction matrices A, which follows e.g. from [SenSl, Definition 2.1, page 46] and [Min 88 , Theorem 
3.2, page 78]. Moreover, Lemma 1.1 stated the existence of a unique positive complex balance equilibrium 
to (2.1) for any given positive initial mass M > 0. Concerning the proof of this part of Lemma 1.1, it 
remains to show the following 


Lemma 2.2 (Unique positive equilibria for strongly connected networks with fixed mass M). 

The first order reaction network JV is strongly connected if and only if the system (2.1) admits a unique 
positive equilibrium for any fixed positive mass M > 0. 


Proof. Sufficiency: Assume that Af is strongly connected. Thanks to the first equivalency in Lemma 
1.1, the reaction matrix A is irreducible. Moreover, for large enough a > 0, we have that A + aE 
is nonnegative in the sense that all of its elements are nonnegative. We can then apply an extension 
of the Perron-Frobenius theorem, see e.g. [SenSl, Theorem 2.6, page 46] or [Per07, Chapter 6.3.1], 
to obtain the existence of a unique positive equilibrium, i.e. a positive right zero-eigenvector X^o = 
(mi.oo, M2,oo, • ■ •, MAT.oo) > 0 Satisfying AXoc = 0 such that ^*.oo = M > 0. 

Necessity: Now assume that (2.1) has a unique positive equilibrium Xoo. Since AX^o = 0 and Xco is 
uniquely determined by the mass conservation, we obtain that dim(ker A) = 1. 

By using a contradiction argument, we assume that J\f is not strongly connected, then the reaction 
matrix A is reducible, i.e. 


A = P^ 




P 


for some permutation matrix P, in which D is irreducible. Choose d to be an eigenvector of D corre¬ 
sponding to zero eigenvalue. Then 








which means that 


positive, P^ 


and Xo 


I is an eigenvector of A corresponding to zero eigenvalue. Since X^o is strictly 
are linear independent, which leads to a contradiction with dim(kerA) = 1. □ 


In the following, we will use the entropy method to study the trend to equilibrium. More precisely, 
for two trajectories X = (ui,U 2 , ■. ■ ,upf) and Y = (vi,V 2 , ■. ■ ,vn) to (2.1), where Y{t) has non-zero 
components for all times t > 0 , we consider the following quadratic relative entropy functional 

£{X\Y){t)=Y, [ ^-^dx. (2.4) 

i=i dn 


The following key Lemma 2.3 provides an explicit expression of the entropy dissipation associated to 
(2.4): 


Lemma 2.3 (Relative entropy dissipation functional). 

Assume that Vi(t) 0 for all i = 1,2,..., N and t > 0. Then, we have 


V{X\Y) 


N 


--£{X\Y) = 2Y,d. 


V v, J 


i N 

dx -1- 





Proof. For convenience we recall that 

N 

dtUi - d^Aui = ^ aijUj 

i=i 


N 

dtVi - diAvi = ^ aijVj, 


and 
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for all I = 1,..., iV. Hence, we compute 


i=l 

^ f 


= / 12—dtUi - ^dtVi dx 

Vi vf J 

I 2^ ^ a,jUj ) - ^ ( d^Av^ + ^ 


N 


Vi 


N 


= E/( 2di—Aui — di—^Avi \ dx + ttijVj j dx 


n\ V, 

(i) 


i=i 


v!\ 

N 




,Ui 


i=i 


(ix 


^_1 ^ Q \ Vi 


N 


* i=i 


A 


= ■ H H ^ 


i=l 


/n 


i=l 


=: Id +Ir. 


(2.5) 


Using integration by parts, we have 


f jil^dx= f ( 2di—Aui — di^Avt^ dx 

Jq Jn V vf J 



Thus, 


Id 




(2.7) 


For the reaction terms Ir, we use an = — X)^i j^i ciji to calculate 




,Ui 


N 


N 


^Vii 
= 2 — 
Vi 


N 

i=ij/* 


dijlij “h diilli I 2 


N 

^ ^ ^ijVj diiVi 


N N ^ 2 / ^ ^ 

“2;~( aijUj — Ui j aijVj — Vi a^, 


i=ij/* 


i=i .i/i 


j=i 


i=io¥* 


A 


— f 2 [aijUj ajiUi) 9 ^jiVi) ) • 

• t ■/■ \ Vi Va 


( 2 . 8 ) 
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Therefore, 


N. . N N , ^ ^2 

^ ^ j Jdx — I ^ ^ ^ ^ I 2 -{^(Xijlij djiUi^ ^jiXi^ j dx 


In 


m , 


i=l j=lj/i 
N 

= E 


(2.9) 


O f \ f \ 

2 {^dijUj CLjiUij 2\^ij'^j ^ji'^i) 

nl '^i 


^UH , N '*^1 / N 

H“ 2 yQijiUi (XijUj J 2 J 

Hj Hj 


dx 


In . 


Ui Uj 

n\Vi Vj 


Vi V. 


2{(Xlj%lj I j {ciijVj i ^ ^ 


N 

= E 

N 

= E 

ui=i;*<j 

^ r / ■ \ ^ 

= “ E / ( S' “ 5" ) 




2{aijUj - ajiUi) - (aijVj - ajiVi) ( — + 


/ J 


dx 


dx 


i,3 = V,i<3 ' 


Vi V, 


By combining (2.5), (2.7) and (2.10), we obtain the result stated in the Lemma. 


( 2 . 10 ) 

□ 


In order to simplify the following calculations, we introduce the difference to the equilibrium 

W := {wi,W 2 , . . .,Wn) = {ui - Ml,oo,W 2 - U 2 ,oo, ....Un- Un,oo) = X - Xoo, 

and remark that thanks to the linearity of the system, the difference W is the solution to (2.1) subject 
to the shifted initial data 


lT(a;, 0) = X(a;, 0) — Xaa , for all a: G Q. 
Note that the total initial mass corresponding to W is zero, i.e. 


i=l 

and that W conserves the zero mass 

N 


AT „ ^ r 

Mw ■= / Wifidx = '^ / (Mi,o(a;) - Ui^oo) da; = 0 , 


2^ / dx = 0, for all t > 0. 


By using the relative entropy dissipation fnnctional derived in Lemma 2.3, we have 


N 


^(fL|-^oo) — 2 ^ ^ j di dx “h ^ ^ (^2j^j,oo “f 


N 


■ 1 Jn lli,oo . . 1 - - - 

1=1 i,j = l;i<j 


Wi Wi 


dx. 


O 111,00 Vj^c 


The following Lemma about entropy-entropy dissipation estimate is the main key to prove the con¬ 
vergence to equilibrium for (2.1). 

Lemma 2.4 (Entropy-Entropy Dissipation Estimate). 

There exists an explicit constant A > 0 depending explicitly on the reaction matrix A, the domain 11, the 
diffusion matrix D and the initial mass M such that 

ViW\Xoo) >\£{W\X^). 

Proof. We divide the proof in several steps: 

Step 1. (Additivity of the relative entropy w.r.t. spatial averages) 

Straightforward calculation leads to 


N 


£{W\X^)=Y^ I ^dx = J2 


N 


\Wi - w^\ 


N 


_ 'Q ^i,oo JQ 

= £(W -W\X^) + £(W\Xoo) 


■dx 


i=l 


\Wi\ 

^i.cn 


( 2 . 11 ) 


where we denote W = {wi, W 2 , •. •, wn) and we recall that Wi = Wi dx for i = 1,..., X due to |11| = 1. 
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Step 2. (Entropy dissipation due to diffusion) 
By using Poincare’s inequality 


we have 


||V/f >Cp|l/-/||2, ior all f € H\n), 
|Vu;,i 2 


( 2 . 12 ) 


N 


N 


-V{W\X^)>Y,d^ / ^^dx>CpY,d. 


-dx 


Q ^i,oo 


i—l ^i,oo 


> Cp min{di,d2, • ■ •, d]si}£{W - W\Xao)- 


(2.13) 


Step 3. (Entropy dissipation due to reactions) 
From (2.11) and (2.13), it remains to control 


£{W\X^)=Y, 


N _2 

Wi 


i=l 


Ui, 


By using Jensen’s inequality we have, recalling that |n| = 1, 

N 


1 


2J(JE|Jtoo) ^ rt ^ ajzlXj oo) [( I 

2 . . 5 -;^/ 7n V^*i,oo Myoo/ 

2 

Wi Wj 


dx 


N 








,00 


dx. 


It then remains to prove that 


1 


N 


rj ^ ^ (,Oiij^j,oo J“ O^ji^i^oo) 




IPi Wj 


Ui. 


00 '-‘' 3,00 


dx > -f 


N _ 2 


i=l 


Ui,, 


(2.14) 


(2.15) 


for some 7 > 0 . Note that if both reactions 5 "^ —> S); and Sj —>■ S'i do not appear in the reaction network. 


then we have a^- = Uji = 0 and thus 


UijUj,Qc> “t“ Uji'lli,Qc — d- 


Hence, the expression 


N 


^ ^ {UijUj,oQ J- CLji'lli,oc) 


Wi Wj 


Ui. 


00 *^ 4,00 


dx 


may not contain all pairs (*, j) with i ^ j. However, the weak reversibility of the network allows to make 
all pairs {i,j) with i ^ j appear in the following sense: There exists an explicit constant ^ > 0 such that 


N 


,00 H” (^ji^i,co^ 


U^ 


,00 '^ 3,00 


N 


dx>^ ^ 


Ui 


,00 '-‘' 3,00 


(2.16) 


Indeed, assume that Uij = Uji = 0 for some i ^ j. Due to the weak reversibility of the network, there 
exists a path from Si to Sj as follows 


Cl _ o °‘3231, c> “3r3r-l_ ^ ^ 

Ji — Jji r Jj2 > ■ ■ ■ > Jj.,. — Oj 

with r > 3 and > 0 for all /c = 2, 3,..., r. Thus, with 

0 <cr= min {ciijUi,oc + 0 ^ 1 %, 00 } < + ^jk-ijk'^jk 

[aij ,aji)ii(0,0);l<i<j<N 2<k<r 


, 00 } 


we have 

r 

^ ^S ^3k3k-i'^jk-i,oo ^jk-ijk'^jk,oo) 


^3k 


Wj 


/c =2 


\ '^jk ,00 

^3k- 

r 

' wj: 

y'^jk.oo 

iM 

II 

a 

> 

( W: 




'31 


^ - 1 \^3u 


00 '^jr ,00 


N — 1 \ui 


2 , 00 '^3,OD 


• (2.17) 


2 
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Since there are less than N(N — l)/2 pairs (*,j) with = Uji = 0, we can repeat this procedure to 
finally get (2.16) with ^ = 2a/{N(N — 1)^). From (2.15) and (2.16), we are left to find a constant 7 > 0 
satisfying 

N 


E 


Wi Wi 


27 ^ Wi 


i-rJ: 


. ^2,00 ^7,00 / 4 -I ^2,00 

with the constraint of the conserved zero total mass 

N 

wf= 0 - 

i=l 

Because of (2.19), 

N N N 

wf' = — y^ wiwj = —2 wiwj. 

i=l ij=l;i^j 

Therefore, we can estimate for C = mini<i<j<jv — —^ 7 :— 


(2.18) 


(2.19) 


( 2 . 20 ) 


N 

E 


Wi Wi 


-*'2,00 ^J:00 


> —2 min ■ 


> min ■ 


1 


i'Cj 


E^^ 

i<j 


oo^j,oo 


Wi Wi 


1 


^ 2,00 ^j,oo 


N 


2<J Ui^QQUj^QQ 


E Wi Wj = min-y^ wi^ > min 

i.^i 7/ - __7/- •jV't 


N _ 2 

Wi 


iKj Ui^ooUj^co 


E 


i<j " i=l 

In conclusion, we have proved (2.18) with 7 = | mini<j ^ 
implies (2.15) and thus completes the proof of this Lemma. 


Ui oo^j,oo ■ - Hz,00 


( 2 . 21 ) 


Ui ooUi c 


-, which in combination with (2.16) 

□ 


Theorem 2.5 (Convergence to Equilibrium). 

Consider (w.l.o.g) a strongly connected subnetwork M of a weakly reversible first order reaction network. 
Assume for Af that the diffusion coefficients di are positive for all i = l.,2,, N, and the initial mass 
M is positive. 

Then, the unique global solution to (2.1) converges to the unique positive equilibrium Xoo in the 
following sense: 

2=1 ^i,oo JQ Ui^ao 

where the constant X > 0 is computed as in Lemma 2.4- 


Proof. From Lemma 2.4 we have 

j^SiX - X^\X^) = -V{X - X^\X,^) <\S{X- X^\X^). 

By Gronwall’s inequality, 

£{X{t) - XooIXoo) < e-^*£iXo - x^\x^), 

and the proof is complete. □ 


Proof of Theorem 1.2. Theorem 1.2 is a direct consequence of Theorem 2.5 and the partition of weakly 
reversible first order reaction network into strongly connected components. □ 

We now turn to the case of degenerate diffusion, where some of the diffusion coefficients di can be 
zero. In the proof of the Theorem 2.5, we have used non-degenerate diffusion of all species in order to 
control distance of the concentrations to their spatial averages (see estimate (2.13)). This procedure 
must thus be adapted in the case of degenerate diffusion. 

It was already proven in [DF07, FLT, MHM] that even if some diffusion coefficients vanish, one can 
still show exponential convergence to equilibrium provided reversible reactions. The technique used in 
these mentioned references is based on the fact that diffusion of one specie, which is connected through 
a reversible reaction with another specie, induces a indirect kind of ” diffusion effect” to the latter specie. 

We will prove that this principle is still valid for weakly reversible reaction networks as considered in 
this section. 
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Theorem 2.6 (Convergence to Equilibrium with Degenerate Diffusion). 

Consider (w.l.o.g) a strongly connected subnetwork M of a weakly reversible first order reaction network. 
Assume for M that the initial mass M is positive. Moreover, assume that at least one diffusion coefficient 
di is positive for some i = 1,2,..., N. 

Then, the solution to (2.1) converges exponentially to equilibrium via the following estimate 


N 

E 


i=l •' ^ i=l 

for some explicit rate A' > 0 which depends explicitly on A, Q, D and M.. 

Proof. We aim for a similar entropy-entropy dissipation inequality as stated in Lemma (2.4), i.e. we 
want to find a constant A' > 0 such that 


^2,00 


-dx < e 


-x't 


N 

E 


\Uifi 


'n 


Ui, 


■ dx. 


V{W\X^) > A' £iW\X,^) = A' [£{W - W|Xoo) + £{W\X^)\. 


( 2 . 22 ) 


Due to the degenerate diffusion, the diffusion part of 'D{W\X^) is insuflicent to control £{W —W\Xoo) as 
in (2.13), since some of diffusion coefficients can be zero. This difficulty can be resolved by quantifying the 
fact that diffusion of one specie is transferred to another species when connected via a weakly reversible 
reaction path. Without loss of generality, we assume that di > 0 and estimate 'D{W\Xoo) by 


D(lT|.X^oo) ^ d\ I dx 4“ ^ ( (aijUj^oQ -t- aj^Ui^oQ^ I ( 

Jn «- 1.00 


By arguments similar to (2.16) and (2.17), we have 


dx. 


N 


’DiWlXoo) > di 


I Cl ^1,00 


—dx + c 


[ 1 

f Wi 


In ' 

V ^2,00 

^j,oo J 


dx. 


(2.23) 


(2.24) 


To control £(W — W\XaQ), we use the following estimate for alH = 2, 3,..., iV: 


f iViPip 
I Ul,oo 


dx ■ 


Wi 


Q \^l,oo Ilz,oo 


dx> (3 


j 

JQ. 


■ dx. 


(2.25) 


with P = ^— min l|: Indeed, thanks to Poincare’s inequality || V/|p > Cp\\f — fW^, we estimate 

for various sufficiently small constants C 


f iViPip 

Iq 


dx 


Wi 


Q \lll,oo '^i,oo 


1 . 

> — mm 
- 2 


= — mm 
2 


= — mm 
2 


Cp 

^ 1,00 

Cp 

^ 1,00 

Cp 

^ 1,00 


dx > 


Wi 


Cp 




|wi - Wl| 
"^ 1,00 
2 

dx 


Wi — Wi ^ Wi 
lll.OO 


Ui, 




dx 


f Cl \^l,oo 111,00 

Wl wi 


Wi 


Wi 


f Cl \lll,oo Hz.oo 111,00 Hi 
2 

Wl 


W, 


> 


Cp 

^l,OD 

Now, thanks to (2.24) and (2.25) 


2u 


1,00 


n \^l,00 ^2,00 


^2 - Wi\ 


1 


dx + — min 


dx 

Cp 

^ 1,00 


n \^ 2,00 ^ 2,00 


■ dx. 


ICl 11^^00 


dx 

(2.26) 


V{W\X^)>mm\^,^\pJ2 


N 


m - Wi\ 


2=1 


Ici 


N 






Wl 


Wi 


Cl \lll,oo 


dx 


> min 


<h C 
N' 2 


P£{W-W\X^) + l Y. 


N 


> min 


'<h e 

iV’ 2 
>\'£{W\X^) 




P£iW -W\X^) + AS{W\X,^) 


f i 

f Wl 


In 

^ “^2,00 

^^‘,00 J 


dx 


(by using (2.18)) 


with A' = min 4 Thus (2.22) is proved and the proof is complete. 


(2.27) 

□ 
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Proof of Theorem 1.3. Theorem 1.3 is a direct consequence of Theorem 2.6 and the partition of weakly 
reversible first order reaction network into strongly connected components. □ 

Remark 2.1. The estimate (2.25) is usually interpreted as follows: the sum of the dissipation due to 
the diffusion of wi and the dissipation caused by the reaction between wi and Wi are bounded below by 
(2.26), which is essentially a diffusion dissipation term of the specie Wi (after having applied Poincare’s 
inequality). In this sense, a ’’diffusion effect” has been transferred onto Wi. 

We remark that while the presented proof for the linear case is straightforward, the proof of an anal¬ 
ogous estimate to (2.25) in nonlinear eases turns out to be quite tricky. Readers are referred to [DF07] 
or [FLT, Lemma 3.6] for more details. 

3. NON-WEAKLY reversible NETWORKS 

In this section, we consider (w.l.o.g.) reaction networks Af which are not weakly reversible, yet form 
one linkage class. Thus, the corresponding directed graph G is connected yet not strongly connected (i.e. 
the underlying undirected graph of G is connected). We will show that in the large time behaviour, each 
specie tends exponentially fast either to zero or to a positive equilibrium value depending on its position 
in the graph representing the network. 

For weakly reversible reaction-diffusion networks (corresponding to strongly connected graphs), it was 
proven in Section 2 that each specie converges exponentially fast to a unique, positive equilibrium value, 
which is given explicitly in terms of the reaction rates and the conserved initial total mass. 

For non weakly reversible reaction networks, however, we will show that while the equilibria are still 
unique and attained exponentially fast, the equilibrium values are in general no longer explicitly given 
but depend on the position in the graph in general and on the history of the concentrations of the 
influencing species in particular. 

Moreover, since non weakly reversible reaction networks (2.1) may no longer have positive equilibria, 
the relative entropy method used in Section 2 is not directly applicable. Nevertheless, we will see that 
the relative entropy and the ideas of the entropy method still play the essential role our analysis of 
non-weakly reversible networks. 

As the large time behaviour of the species depend on their position within the network, we need to 
first state some important properties of the graph G. The following Lemmas 3.1 and 3.2 are well known 
in graph theory. We refer the reader to the book [BJG08] for a reference. 

Lemma 3.1 (Strongly connected components form acyclic graphs G*^). 

Let G be a directed graph which is connected, that is the underlying undirected graph of G is connected, 
but not strongly connected such that the graph G contains at least r > 2 strongly connect components, 
which we shall denote 6 j/ Gi, G 2 ,..., G^. Thus, we can define a directed graph G^ of strongly connected 
components as follows 

- G^ has as nodes the r strongly connected components Gi, G 2 ,..., G^, 

- for two nodes Gi and Gj of G^, if there exists a reaction Gi 9 Sk Si € Gj with aik > 0, 
then we define a directed edge Gi —>■ Gj of G^. 

Then, the directed graph G*^ is acyclic, that is G^ does not contain any cycles. 

Proof. The proof can be found in e.g. [BJG08, Ghapter 1] and shows that if G*^ would contain a cycle 
then this cycle should have been contained in a strongly connected component in the first place. □ 

Lemma 3.2 (Topological order of acyclic graphs, [BJG08, Chapter 1]). 

There exists a reordering of the nodes of G^ in such a way that for all direct edges Gi —>■ Gj we always 
have i < j. 

From now on, we will always consider topologically ordered graphs G^. For each i = 1,2,..., iV, we 
denote by Ni the number of species belonging to Gi. For notational convenience later on, we shall set 
L[0] = 0 and introduce the cumulative number L[i] of the species contained in all strongly connected 
components up to Gi, i.e. 

L\i] = Ni-^ N 2 +...-I Ni for alH = 1, 2,..., r. (3.1) 

We then reorder the species of the network Af in such a order that the species belong to the component 
Gi are S'L[i_i]-Hi, S'L[i_i]-H 2 , ■ • ■, for alH = 1,2,..., iV. 

Each component Gi belongs to one of the following three types: 


ENTROPY METHOD FOR REACTION-DIFFUSION SYSTEMS 


18 


• Source component: Ci is a source component if there is no in-flow to Cj, i.e. there does not exist 
an edge Ci ^ Sk ^ Sj G Ci, 

• Target component: Ci is a target component if there is no out-flow from Ci, i.e. there does not 
exist an edge Ci B Sk ^ Sj ^ Ci, 

• Transmission component: If Ci is neither a source component nor a target component, then Ci 
is called a transmission component. 


The above classification of strongly connected components greatly improves the notation of the cor¬ 
responding dynamics, which quantifies the behaviour of the species belonging to the three types of 
components. In the following, we denote by Xi = (ui[i_i]+i, Ui[i_i]+ 2 ,..., the concentration 

vector of the species belonging to Ci. 

The evolution of the species belonging to a component Ci depends on the type of Cp. 

(i) For a source component Ci, the system for Xi is of the form 

r dtX, - D,AX, = A,X, - Fr*X„ xGil, t>0, 

} d^x, = 0, xGdn, t>0, (3.2) 

0 ) = x G fl, 

where the diffusion matrix Di is 

A = diag(dL[,_i]+i,dL[,_i]+2,... ,dL[i]) G (3.3) 

the reaction matrix Ai is 

Ai = (ai[i_i]+fc_i[j_i]+^)i</c,£<iVi G (3.4) 

and the out flow matrix is defined as 


Fr* = diag(A[,_i]+i, /i[,_i]+2 ,..., Ah) G (3.5) 

with 

N 

f L[i — l]+k ^ ^ • 5 

t=L[i] + l 

where the lower summation index L[i\ + 1 follows for the topological order of the graph C'^. 

Roughly speaking, A[i_i]_|_fc is the sum of all the reaction rates from the specie Afi-ij-i-fc to 
species outside of Ci. It may happen that A[i_i]_|_fc = 0 for some k = 1,2,..., N, but there 
exists at least one ko such that fLli-i]+ko > 0 since Ci is a source component. 


(ii) If Ci is a transmission component, the system for Xi writes as 


idtX,-D,AX,=Fl^+A,Xi-F°'^^X„ x G n, t > 0, 

< d,^Xi =0, X G dll, t > 0, (3.6) 

[ A(a:, 0) = A,o(a;), x G Tt, 


where the diffusion matrix Di, the reaction matrix Ai and the out flow matrixF°“* are defined 
as above in (3.3), (3.4) and (3.5), respectively. The in-flow vector F™ is defined by 


A' 


^ L [ i - l ]-|-2 

V zm ) 


L[i-1] 

with ^ ■ 

k=l 


(3.7) 


We remark that by studying all components Ci within the topological order of , the dynamics 
of the previous components Ci,C 2 ,. ■., A-i is already known at the time we analyse the compo¬ 
nent Ci. Thus, in system (3.6) the in-flow vector Jy" can be considered as a given external in-flow. 


(hi) If Ci is a target component, we can write 

r dtX, - AAA = + A,X„ xGlt, t>0, 

< di,Xi = 0, X G dll, t > 0, (3.8) 

[ A(a;,0) = W.o(a;), x G 11, 

where the reaction matrix Ai and the in-flow Jq*" defined in the same way as above in (3.4) 
and (3.7). 
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By modifying the relative entropy method in Section 2, we obtain the 


Proof of Theorem 1.4-- Since the ongoing outflow vanishes the mass of all source components and sub¬ 
sequently all transmission components, the corresponding equilibrium values are expected to be zero 
and the relative entropy method used for weakly reversible networks is not directly applicable here. We 
instead introduce a concept of ” artificial equilibrium states with normalised mass” for these components, 
which allows to derive a quadratic entropy-like functional, which can be proved to decay exponentially. 
Due to their different dynamics, we have to distinguish the two cases: Ci is a source component and Ci 
is a transmission component. 

The aim of the proof is to show that if Ci is a source or a transmission component then for all 
k=\,...,Ni, 

||■UL[i_l]+fe(<)||| 2 (^) < K,e~^'\ for all t > 0, (3.9) 

for explicit constants Ki > 0 and > 0. 

In order to simplify the notation, we shall denote 

Vk = UL[,_i]+k, and bk,e = aL[i_i]+k,Lli-i]+i, for all 1 < k,£ < Ni. (3.10) 

Then, the concentration vector Xi and the reaction matrix Ai can be rewritten as 
Xi = {vi,V 2 ,...,VNi) and A^ = {bk,e)i<k,i<Ni- 
Note that the index i for the component Ci is fixed. 


Case 1: Ci is a source component. 

We recall the corresponding system from (3.2) 


dtX, - D,xx, = AX, - Fr*x„ xen, t>o, 

d„Xi = 0, X € dfl, t > 0, 

X,{x,0) = xGfl. 


(3.11) 


We now introduce an artificial equilibrium state Xi^^o = ('Ci,oo, 112,00 ■ ■ ■ ,VNi,oo)'^ with normalised mass 
to (3.11), which is defined as the solution of the system 


I AiXi^oc — 0, 

I 111, 00 T P2,oo -f . . ■ -f VNi^oo — 1- 


(3.12) 


It follows from Lemma 2.2 that there exists a unique positive solution W,oo to (3.12). Here we notice 
that Xi^oo balances all reactions within Ci while the total mass contained in Xi^^a is normalised to one. 

In the following we will study the evolution of the quadratic entropy-like functional 

Ni 


£iX,\X,,oo) = 


\Vk\ 


■dx. 




f2 '^k,oc 


(3.13) 


By similar calculations as in Lemma 2.3, we obtain the time derivative of this quadratic functional 


dt 


Ni 


= 2j2d 




k^l 


f iVPfcP 

O '^k,oc 


■dx 


Ni f f V 

+ V' (bk,ive ,00 + bi^kVk , 00 ) / ( I 

Jn\Vk,oc J 


dx 


k,^—l-,k<i 
Ni 

' 2 ^ ^ fLli-l]+k 

fe=l 


\m 


-dx. 


n ^fc,oo 


(3.14) 
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We remark that since Ci is a source component, there exists an index fcp S {1, 2,..., Ni] such that the 
out-flow fL[i-i]+ka > 0 is positive. Then, an estimate similar to (2.16) gives for various constants C 


Ni 


P(W|^z,oo)>e 

k,£—l\k<i 


/ f """ 

Vi 

Y 

Jq ^ '^k,OQ 

'^1,00 

) 


, Of f 

1 “1“ L[i—l]-\-kQ I dx 

,oo / Jq '^ko,oo 


> min 


inU/2,/i[,_i]+fc„/2iV,} E / (^ 

I JL[i — l]-\-ko I 

JCl '^ko,oo 

p \ 12 

>\y = A,£:(X,|X,.oo) 

Jn ve,oG 


Vi Vka 


oo VkQjOO 


-s 


+ 


\vko 


Vko 


dx 


(3.15) 


with Xi = min{^/4,It follows that 


-f(X,|X,,oo) = -P(X,|W,oo) < -\i£{Xi\Xi^i^), 

dt 


and thus 


Ni 

E 


k/=( 0 l" 


dx = £iXi{t)\Xi,oo) < e-^^*f(W.o|^^,oo) 


or equivalently 

\\uL<i_i]+kitW < e~^'*£{Xi^o\Xi^i^) max {ui^oo} for alH > 0 , for all fc = 1 , 2 ,..., W, 

l<i<Ni 

which proves (3.9) with Ki = £{Xifi\Xi^ao) inaxi<i<jv^{ui_oo} in the case Ci is a source component. 
Case 2: Ci is a transmission component. 

By recalling that the components Ci are topologically ordered, we can assume without loss of generality 
that Ui, with £ = 1 , 2 ,..., L[z — 1 ], obeys the following exponential decay 


\uiit)r < K*e 


* ^—X*t 


'= 1,2,L[i — 1], for all t > 0. 


for 0 < A* = min Afe and K* = max Ki. We also recall the system for Ci, 

l<k<i-l l<k<i-l 

' dtXi - DiAXi = J-™ + AXi - Fr*Xi, xeA t>o, 
dijXi = 0, X € dH, t > 0, 

Xi{x,0) = Xi^o{x), X G ft, 


(3.16) 


(3.17) 


with J"*” is defined as (3.7). Denote by Xi^oo = ('Pi.oo, ■ ■ ■,'PAi.oo)^ the artificial equilibrium state of 
(3.17), which is the unique positive solution to 


J oo — 0 , 

("^1,00 4“ V2^oo 4“ ... “t“ rJ7Vi,oo — !■ 

Again, we can compute the time derivative of 


Ni 


£(x,iw,oo) = 51 


Vk 


■dx 


as 


P(W|^*.oo) =-:w£(X*,^z.oo) 
dt 


— 2 5 ^ C^L[i- 


f |V^^fcP 

O '^k,oo 


fc=l ' 


Ni f f V} 

dx “h ^ ^ (^/c,£'^£,oo H“ ^£,fc^/c,oo) / ( 

k,e=l-,k<i 

+ 25^/L[i_i]_^fc / -^dx-25^ / ^ aL[i-i]+k,eUi]dx. 

k=l JnVk,oo ^^^Jn\Vk,ca / 


n Vk,i 


(3.18) 


(3.19) 


Vk Vi 


00 '^£,00 


dx 


(3.20) 
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Because Ci is a transmission component, there exists an index fcg S Ni} such that fLli-i]+ko > 0 

is positive. In comparison to (3.14), the dissipation ViXilXi^oo) in (3.20) has the additional term 


Ni r / ^ ^F- 1 ] \ 


to be estimated. Thanks to the decay (3.16) of U(, we can estimate 

Ni „ / 


iVi „ y -^J \ 


Ni T[i-1] 

<2^ ^ aL[i_i]+k,e 


k=l l=\ 


f 

Vk 

/ 

Ul 

/a 

^fc,oo 


dx 


Nj „ I 12 

</L[i_l]+feo^ / + K ^ ||uff 

p \ 12 

< ^ (3.21) 

./n ^fc,oo 

with K = NiL[i — 1] max{a^^}/(/i[i_i]+fe„ min{ufc_oo})- Then, with the help of (3.21), we estimate 


V{Xi\Xi^i,o) > ^ ibk,m^ac + be^kVk,oo) [ (— - dx + fL[i-i]+ko [ ^—^^dx - kK* 

k^e^k<e Jn\Vk,oo vi,ooJ JnVko.oo 


,-A*t 


and similarly to (3.15), we obtain for A = min{^/4, 

P(X,|X,,oo) > Af(X,|X,,oo) - KK*e-^^\ (3.22) 

From (3.22), we can use the classic Gronwall lemma to have 

E{Xi{t)\Xi^i^) < 

with \i = min{A,A*} and Ki = 2 max{£i(Afi^o|-^i,oo), nIF*}, which ends the proof in the case that Ci is 
a transmission component. □ 


For a target component, we need to define its corresponding equilibrium state. This equilibrium state 
balances the reactions within the component and has as total mass the sum of the initial total mass of 
the target component plus the total ’’injected mass” from the other components. In general, the injected 
mass will not be given explicitly but depend on the time evolution of the influences species prior to Ci 
in terms of the topological order. 


Lemma 3.3 (Equilibrium state of target components). 

For each target component Ci, if 

Ni Ni L[i-1] 

aL[i-i]+k,t / ut{s)ds > 0 (3.23) 

fe=l fe=l ^=1 g 

holds, then there exists a unique positive equilibrium state Xi^ao = ("Ci.oo)'P 2 ,oo 5 • ■ • ,VNi,oo) satisfying 


— O 5 

Ni Ni Ni - 1-00 (3.24) 

IIfc,oo — l]-t-fc,0 4“ g^ ^L[i—l\ + k,t f U^(s)ds. 

k=l k=l k=l e=l 0 

Otherwise, if the sum (3.23) should be zero, then the initial and the total injected mass into the target 
component Ci is zero and the concentrations of the target component Ci remain zero of all times. 


Proof. By (3.16) we have for all £ = 1,2,..., L[i — 1] that ||u^(t)p < K*e ^ *. Thus, Jensen’s inequality 
yields 


^ + 00 


^ + 00 


Ui{s)ds < I ||w(s)|li 2 ^(o)ds < K* 


n+oc 


e-^^ds = 


2K* 


(3.25) 


and the right hand side of the second equation in (3.24) is finite. Therefore, the existence of a unique 
Xi^ao satisfying (3.24) follows from Lemma 2.2. □ 
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Remark 3.1. The positive sign in assumption (3.23) ensures that either initially or during the ongoing 
reactions positive mass is present/injected into the component Ci. When this assumption does not hold, 
then the target component does not possess a positive equilibrium and all of its concentrations remain 
zero for all times. For example, consider the network 

GSl 

-^^3 

0-21 

S2 -. ^4 

£*42 

when the initial data of all species are zero except S 3 . In this case, the target component { 5 ' 4 } will not 
ever receive any mass, and thus remains zero for all t > 0 . 

We now begin the 


Proof of Theorem 1.5. With the notations introduced in (3.1) and Lemma 3.3, we identify the indexes in 
the statement of Theorem 1.5 as ik = L[i — V\ + k and the equilibrium state 11.4,00 = Vk,oo for fc = 1,..., Ni. 
The aim now is to prove for all A: = 1,..., W, 

\\vk{t) - Ufe,oo||i 2 (o) < Kie~^** for alH > 0 

for some explicit constants Ki > 0 and Ai > 0. 


We recall the system for a target component Ci, 


where 


dtXi - DiAX, 

_ 

+ AiXi, 

X G £ 1 , 

t > 0, 

d,X, = 0, 




X G d£l, 

t > 0, 

7 ^ 

0 

II 


(a:), 


X G £ 1 , 






L[i-1] 


ZL[i-l]+2 


with 


= E 

—l]+£,fc 

V zm / 




k^l 



(3.26) 


Note that the total mass of Ci is not conserved but increases in time due to the in-flow vector Jq™. 
To compute the total mass of Ci at a time t > 0, we sum up all the equations of (3.26) then integrating 
over n, 


fL 

dt 


Ni Ni Ni 

Y,UL[i-l] + k{t) = ^ZL[,_i] + fc(t) = ^ ^ ^L[i — l]+k,t '^i{t') 
k=l k=l k=l i=l 


thanks to the homogeneous Neumann boundary condition and the fact that (1,..., 1)^ is a left eigenvector 
with eigenvalue zero of Ai since Ai is a reaction matrix. Thus, we have 


Ni Ni Ni L[i-l] . 

^?lL[,_l]+fe(t) = ^UL[i_i] + fc,o + X! / ue{s)ds. (3.27) 

fc=i fe=i k=i e=i 


Given that the right hand side of (3.27) should be zero for all times t > 0, then ?Ij;,[j_i]_|_fe(t) = 0 for all 
k = 1,..., Ni and for alH > 0 and Xi^oo = 0 and the statement of the Theorem holds trivially. 

Otherwise, if the right hand side of (3.27) is positive for some time t > 0, then assumption (3.23) is 
satisfied an Xi^^o is a positive equilibrium. Recalling the change of notation Vk = in (3.10), we 

denote by 

'^kit^ — '^k{t') Vji^co — ^L[I—l]-t-fc (a) 

the distance from to its corresponding equilibrium state for all k = 1 , 2 ,...,W- It implies 

that {wk)k=i,...,Ni solves the system (3.26) subject to the initial data Wkp = UL[i_i]-|-fc,o ~ 'Pfc.oo for ah 
k = 1,2,..., Ni. We define Wi = (wi, W 2 ,..., wjVi) and consider the relative entropy-like functional 


Ni 




{Wkl 


Ni 


JQ, '^k,oo 


■dx = J2 


= l * 


\Wk - Wk\ 
'^k.oo 


-dx + 


Ni _ 2 

Wk 




—: £i + £ 2 . 


(3.28) 
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By using again arguments of Lemma 2.3, we calculate the entropy dissipation 


'D{W,\X,^oo) = -J^S{W^\X,^oo) 


Ni 






2=1 


f \^Wk\- 

O '^k,oo 


f f w 

dx “h ^ (bk,£'^i,(yo “t“ ^^,/c'^/c,cso) / ( 

k,e^k<e 


Wk Wj 

OD “^£,00 


Ni 1] 

- 2 ^ ^ ®T[i-l]+fc,^ 
fc=i 

For the last term of (3.29), we estimate 

Ni L[i-1] 


Wk 


Q ^fc,oo 


Ui dx 


dx 


(3.29) 


2EE 0‘Lli-l]+k,e 
k=l t=l 

Ni L[i-i\ 

< 2 ^ ^ ak[i-i]+k,t. 


Wk 


Q ^fc,oo 


ui dx 


fe=l £=1 
Ni 


< Cp E ^L[2-l] + fc 

fc=l 

Ni 

— ^ ^ *^L[ 2 —l] + fc 


/ \Wk - Wk\ 
Iq '^k,oo 

■|2 


Ni L'l'i 1] |___—I 

\ug\dx + 2 EE ^L[z— 


fc=l {=1 


'^k,oo 


__ L[i-l\ Ni - 2 L[i-1] 

Wk — Wk\ , II ||2 Wk -2 

dx + Ni 2_^ +K2 2^-h K3 2.^ Ui 

, T '^k.oo 

k—1 ’ 


jQ ‘^fc,oo 

iVwfcp 


fci 




fe=i 


Q ^fc,oo 


iVi _ 2 

dx + K 2 ^ -1- (ki + K'i)K*e 


* _ — A*i 


fe=l 




(3.30) 


with 


— 1] max{af ■} NiL[i — 1] max{af } 

U 1 = -N, -r-j-^-y, K 2 = max{ufc,oo}, ^3 = -, 

Op min-lapu—ii-^fc'L’fc^ooj ^ ^2^fc,oo 

k 

where /^2 is chosen in such a way that the last step of the below estimate (3.35) is fulfilled, and we have 
used ||u^(t)|p < K*e~^ * for all £ = 1,..., L[z — 1] in the last estimate. By inserting (3.30) into (3.29), 
we obtain 


Ni 


^?(VF^|X,.oo) > E [i —1]+A; 


f |Vw/e|^ 

Q 'Cfc^oo 


Ni 


■dx + E {bk,iVi 

,oo + bi^kVk ,oo ) 


k,£—l\k<l 


Wk Wi 


Vk,i 


Vi, 


Ni _2 

Wk 


- N2 y —^ - (ki + K3)K*e-^^* =: Pi + Pa 

fc=i 


(3.31) 


where Pi is the term containing the gradients and Pa is the rest of the right hand side. It follows from 
Poincare’s inequality that 

Ni „ 1 , Ni 


T^l > E ^L[i- 




f IViPfch 

n '^k,oo 


dx > Cp'^dL[i_ 


i]-i-fe 


/ \Wk - Wk\ 
Q Vk^oo 


■dx > K 4 £i 


(3.32) 


with K 4 = Cp min{di^ri_i]_|_fc}. To control £ 2 , we use arguments similar to Step 3 in the proof of Lemma 

k 

2 . 4 . First, by using (3.27), we have the total mass of {wk)i<k<Ni is computed as, 

Ni Ni Ni Ni Ni L[i-1] . 

E= E mi-i]+kit) - E = E mi-i]+k,o + E E ^L[i—i]-\-k,i I uii^s^ds 

k—\ k—\ k—\ k—\ k—"]. £—\ 

Ni Ni L[i-l] , ^ 

-E“i[*-ii+M-E E UL[i — l]+k,l I Ui(^s')ds 

k=l k=l 1=1 ■'^ 

Ni T[i- 1 ] 


Hence, 


Ni Ni 

y wfe w=- E 

k,i—l-,k<£ k,£—l\k^£ 


l\i /*+oo 

= -E E UL[i-i]+k,e / ue{s)ds =: (3.33) 

k=i e=i 

Ni Ni Ni 

WkWi = '^Wk'^ - y WkWi = '^Wk'^ - 


k=l 


k,e=i 


k=l 


-2 


(3.34) 
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By using (2.16) and (3.34), we estimate 


Ni 




Wk 


Wl 




Ni 


- At2 E — - -+ K3)K*e ^ * 

'^kjOO 




Ni 


Ni _ 2 

Wk 


> -2^ max{z;fc ooV^.oo} E Wfc - ^2 E — - -(^i + i^3)K*e 

k<i ^^ ^^ - - 




' Ni 


k=l 


Vk. 


Ni _ 2 

Wk 


= ^ Taa,y.{vk,ocV(, ao} E - ^2 E-(^i + i^3)K*e 

k<i \ I '?’»_ 




Ni 


\k^l 


fe=l 


'^k,oo 


> max{iifc oo} E — - -maxjufe oo^'^,oo}<5^ - (^i + K3)K*e 

2 k * ^ - - h<'P 


-x*t 


/c=l 


Vk 


k<l 


for e > 0 is sufficiently small. It follows from (3.33) and ui < ||m^|| < \fK*e ^ that 

Lli-'i] / ^+oo X 2 


-^ 1 * / p-\-CiO 

< iViL[z — 1] inaxla^ } E f / W{s)di 


s < K^e 


— X*t 


with K 4 = 4:K*NiL[i — 1]^max{af }(A*) Hence, (3.35) implies that 

i<j 


(3.35) 


'^2 > max{ufe,oo} E - - max{K;4^ maxjufe oo'Cf.oo}, (ki + K 3 )K*}e ^ * = ^,362 - n^e ^ *. (3.36) 

2 k ^ Vk,oo k<i 

k—1 ’ 

Combining (3.36) and (3.32) yields 

V{W,\X,^^) > min{K4,/C5K(W',|X,,oo) - (3.37) 

Therefore, by applying a classic Gronwall lemma, 

£{W^{t)\Xi^^) < forallt>0 (3.38) 

with \i = min{^^ 4 , A*} and Ki = 2 max{^(Ar^^o|^z,oo), /^e}- This completes the proof of the Theorem. 

□ 
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